Neutrino Emission from Gamma-Ray Burst Fireballs, Revised 
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We review the neutrino flux from gamma-ray bursts, which is estimated from gamma-ray obser- 
vations and used for the interpretation of recent IceCube data, from a particle physics perspective. 
We numerically calculate the neutrino flux for the same astrophysical assumptions as the analytical 
flreball neutrino model, including the dominant pion and kaon production modes, flavor mixing, 
and magnetic field effects on the secondary muons, pions, and kaons. We demonstrate that taking 
into account the full energy dependencies of all spectra, the normalization of the expected neutrino 
flux reduces by about one order of magnitude and the spectrum shifts to higher energies, where we 
can pin down the exact origin of the discrepancies by the re-computation of the analytical models. 
We also reproduce the IceCube-40 analysis for exactly the same bursts and same assumptions and 
illustrate the impact of uncertainties. We conclude that the baryonic loading of the flreballs, which 
is an important control parameter for the emission of cosmic rays, can be constrained signiflcantly 
with the full-scale experiment after about ten years. 



If gamma-ray bursts (GRBs) are sources of ultra-high 
energy cosmic rays (UHECRs), they should also lead to 
neutrino production [1] [2]. The IceCube [3] (IC) neu- 
trino telescope has, for the first time, significantly con- 
strained the neutrino flux to below the expectations from 
gamma-ray and cosmic ray observations [U E] ; see also 
Ref. [6] for a fit to cosmic ray data. In particular, the 
method to compute the expected neutrino flux from the 
gamma-ray fluence in the internal shock model has been 
derived in Refs. [TJ H]; we refer by "IC Fireball (neu- 
trino) Calculation" (IC-FC) exactly to the IceCube ver- 
sion in Ref. [8| , which is based on Ref. [7] . On the other 
hand, it has been clear from numerical calculations that 
there are limitations to the analytical method from the 
particle physics perspective, see, for instance, Refs. 
I15j for the impact of additional pion/kaon production 
modes, flavor mixing, and magnetic field effects on the 
secondary muons, pions, and kaons. For example, nor- 
malizing the proton and photon densities in the source 
to the Waxman-Bahcall (WB) ORB flux |2], it has been 
demonstrated in Ref. [13] that the combination of these 
effects modifies the shape significantly, and increases the 
normalization by a factor of three to four. So obviously 
there has been increasing tension between theory and ob- 
servation, which has challenged the paradigm that GRBs 
are the sources of the UHECR. We study the connection 
between gamma-ray observations and neutrinos by re- 
interpreting the IC40 data with a numerical model based 
on exactly the same assumptions, same parameters, and 
same bursts, i.e., without changing the astrophysical in- 
gredients. However, we include the additional multi-pion, 
kaon, and neutron production modes, the synchrotron 
losses of the secondaries, adiabatic cooling, and the full 
energy dependence of the spectra. Compared to Ref. [M] , 
we do not normalize the neutrino flux to the WB flux, 
but to the actually observed photon fluence. That is, 
the photon density in the source is obtained from E}^° 
following the gamma-ray observation, the magnetic field 



is obtained from energy equipartition between electrons 
and magnetic field, and the proton density is assumed 
to follow an injection spectrum with the normal- 
ization determined by the baryonic loading. Since we 
find significant discrepancies in normalization and shape 
compared to the analytical models, we pin down the dif- 
ferences by re-computations of the (original) analytical 
models in Refs. [T1[71[S]. 

First of all, consider the analytical method IC-FC in 
Ref. [S] (see App. A therein), used for the IceCube anal- 
yses. At the source, protons, injected with an E^^ spec- 
trum, are assumed to collide with target photons with 
a broken power law spectrum which comes from the 
gamma-ray observation on a burst-by-burst basis. The 
neutrino spectrum is assumed to have two breaks, one 
from the photon spectrum, and one from the cooling of 
the secondaries. The normalization of the neutrino flu- 
ence is computed from the photon fluence with 
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where {xp^^y) ~ 0.2 is the (average) fraction of proton 
energy going into a pion per interaction, /e is the frac- 
tion of the total energy in electrons compared to the to- 
tal energy in protons {^/fe- baryonic loading), /^r is the 
pion production efficiency, Xpj = l/{n^(J\) is the pro- 
ton mean free path, and Ai? is the shell width. On the 
other hand, our Numerical Fireball Calculation (NFC) is 
described in detail in Ref. [T5] (model "FB-D" therein). 
In short, the model relies on the proton and photon den- 
sities within the source. Once these spectra are fixed, 
the rest is just particle physics, where the effect of syn- 
chrotron cooling, adiabatic cooling, and decay of the sec- 
ondaries is explicitely computed (compared to the analyt- 
ical approach) . The normalization of the photon density 
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FIG. 1: Left panel: Shape and normalization modification for the muon neutrinos from GRB080603A between IC-FC and 
the Revised (analytical) Fireball Calculation (RFC), where the individual contributions are described in the main text (burst 
parameters: photon spectrum with a = 1, /3 = 2, ejj, = 0.2 MeV, fluence = 5.5 • 10~^ erg cm"'^, F = 10'^'^, = 0.01s, 
Tgo = 180 s, z — 1.69, = 6.0 • 10^^ erg s~^ |16| : B ~ 90 kC, as a consequence of the assumed energy equipartition 

fractions and parameters; see Eq. (17) of Ref. [H]). Right panel: Numerical reproduction of RFC using the WB A-resonance 
approximation [T], and the finally obtained calculation NFC including the additional pion and kaon production modes. 



is obtained by calculating the equivalent energy from the 
measured photon fluence 

^T-T^ / de,e,F,{e,). (2) 
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Then the proton and magnetic (energy) density normal- 
izations are obtained by the usual energy partition as- 
sumptions with the same parameters as in Refs. [H [8], 
and the same assumptions for the geometry of the fire- 
ball. Note that the IceCube analysis is based on a number 
of bursts for which the neutrino flux is "stacked" , since 
the expected neutrino signal from one burst is too small. 
Therefore, this computation has to be performed for each 
burst individually. 

Let us now compare the results of IC-FC and NFC 
by using a simplified analytical version of the numerical 
code, based on the photo-meson production in Ref. [T] 
("WB A-approx"), and by re-computing the analytical 
models. The main difference has been identified to be 
spectral effects: while the analytical computations often 
rely on estimates using a particular energy (e.g., the pho- 
ton break energy) , the numerical code takes into account 
the full energy dependencies automatically. In Fig. [l] 
we illustrate this in the left panel for one specific exam- 
ple from the IC40 analysis producing a result similar to 
Refs. PQH]. The curve IC-FC shows the analytical ex- 
pectation for the chosen parameter set. As a first step, 
the shape is revised (curve "shape revised") by includ- 
ing a shift of the first break (correction of threshold of 
photohadronic interactions in Eq. (A3) of Ref. j^, see 
Ref. [13], or missing factor in Eq. (3) of Ref. [T]), the fact 
that there are two different cooling breaks for muons and 
pions (including flavor mixing), and a factor of 1 -I- z 



from the effect of the cosmic expansion on the variability 
timescale. As the next step, the correction c/^ to the 
pion production efficiency contains: fc^ (energy of all 
photons approximated by break energy, whereas photons 
distributed according to the photon spectrum; coming 
from Eq. (A13) in [7]), /~ — 0.69 (rounding error in Eq. 
(A15) in [7J), and ~ 2/3 (from neglecting the width of 
the A-resonance in Xp^ instead of using the interaction 
rate; after Eq. (A12) in [7|, but included in Eq. (3) of 
Ref. [T]). The factor cs corrects for energy losses of the 
secondaries and the energy-dependence of the mean free 
path of the protons, see Eq. (11) in Ref. [17] and discus- 
sion therein. Note that this factor is somewhat model- 
dependent because the energy in protons, computed from 
the energy partition and baryonic loading, depends loga- 
rithmically on the minimal and maximal proton energies 
(for an injection spectrum), whereas in Eq. ([l]) only 
the part relevant for neutrino production is taken into ac- 
count. To illustrate that, the dotted curve in Fig. [TJ left, 
shows the result in the extreme case that the minimal 
proton energy coincides with the photo-meson produc- 
tion threshold. Note that fc^ and cs strongly depend 
on the photon spectral indices, and vary from burst to 
burst. Surprisingly, all the corrections go into the same 
direction, which means that the approximations in the 
analytical model were probably a bit on the optimistic 
side. The result can be regarded as Revised Fireball (neu- 
trino) Calculation (RFC). In the right panel of Fig. [T] 
it is shown that RFC matches the numerical result for 
the same assumptions for the photohadronic interactions 
("WB A-approx", as in Ref. [I]) very well. Taking into 
account the additional multi-pion and kaon production 
modes, similar to Ref. [M], the flux increases again, and 
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FIG. 2: Reproduction of the IC-FC prediction for the neu- 
trino (differential) fluence E^F^, compared to the correspond- 
ing IC40 hmit (hght/blue curves; 90% CL). In addition, our 
numerical prediction NFC is compared to the corresponding 
IceCube limit for exactly the same bursts and assumptions 
(black curves). Compare to Fig. 2 in Ref. [4]. 

the final (numerical) result NFC is obtained. In this case, 
the normalization deviates about one order of magnitude 
from the analytical prediction IC-FC, and the shape is 
significantly different, shifted to higher energies. Note 
that we have chosen one analytical method IC-FC for 
the comparison, whereas the detailed comparison to an- 
other method, such as Ref. [T], will depend on the specific 
approximations of the analytical method (whereas NFC 
does not depend on these). 

As the next step, we reproduce the IC40 analysis from 
Ref. [1], based on 117 bursts, using the same neutrino ef- 
fective area and same assumptions, bursts, and parame- 
ters [12]. The result is shown in Fig. [2] (light /blue curves), 
where the dashed curve shows the IC-FC prediction for 
the neutrino flux and the solid curve the corresponding 
IC40 limit. In this case, the bound is below the predic- 
tion, and the original model is under tension. Our result 
is shown as black curves: the prediction is about one or- 
der of magnitude below the limit corresponding to this 
flux shape. This qualitatively different result means that 
IceCube has not yet reached the level where it tests the 
parameters chosen for the fireball model. 

In order to obtain conclusions on the cosmic-ray con- 
nection, or to compare the results from different experi- 
ments, the extrapolation of the fluence to a quasi-diffuse 
flux is needed. It depends on the number of bursts ex- 
pected per year, where 667 has been used [3]. We show in 
Fig. [3] our quasi-diffuse flux prediction ("GRB, all") to- 
gether with the IC40 Hmit, the combined IC59-I-40 limit 
(which has a different flux shape), and an extrapolated 
IC86 limit. In addition, we show different regions and 
curves to illustrate the size of several model- or method- 



specific additional "systematical errors": the statistical 
error coming from the extrapolation from a few bursts 
to the quasi-diffuse flux (for 117 bursts, estimated and 
obtained from Ref. [TS]) and the "astrophysical uncer- 
tainty" for this particular model (envelope of the follow- 
ing independent variations around the assumptions for 
the IceCube analysis: variability timescale by one or- 
der of magnitude [0.001s . . . 0.1s for long bursts], F be- 
tween 200 and 500, proton injection index between 1.8 
and 2.2, and ee/cs, energy in electrons versus magnetic 
field, between 0.1 and 10). As one can read off from this 
figure, neither IC40 nor IC59-I-40 can reach the predicted 
fluxes, even in the most optimistic cases; compared to 
IC59-I-40, a factor of two higher statistics is needed to 
reach the nominal prediction. However, the full scale 
IceCube experiment, operated over about 10 years (ex- 
trapolation), will finally find the GRB neutrinos or sig- 
nificantly constrain the model unless, for instance, the 
number ratio between F > 500 and F ^ 300 bursts (or 
corresponding collision radii) is larger than seven for fixed 
t^, as it can be easily shown. Note that our given as- 
trophysical uncertainty is less model-dependent than the 
one in Ref. , since it does not rely on the origin of the 
target photons, but it includes the effects of synchrotron 
losses. 

We have deliberately omitted one variable from this 
discussion: the baryonic loading l//e, which directly re- 
scales the neutrino flux prediction, as illustrated by the 
arrow in Fig. [3] and as it can be read off from Eq. ([I]) . 
The choice of this parameter is often consistent with a 
coherent picture among cosmic ray, gamma-ray, and neu- 
trino fluxes if the GRBs are the sources of the UHECR, 
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FIG. 3: Prediction of the quasi-diffuse flux (NFC), including 
the estimates for several model- or method-specific system- 
atical uncertainties (see main text). In addition, the IC40 
limit is shown, and two expectations are shown for com- 
parison (IC59-I-40 from Ref. [5] and IC86 extrapolated for 
Al^^'^ ~ 3 X A^g*" from IC40; see, e.g., Ref. [H]). 
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i.e., treating it at the same level as the other parameters 
would change the logic of the astrophysical picture. For 
example, if the protons are magnetically confined within 
the source, the cosmic ray flux will be proportional to 
baryonic loading and the photohadronic interaction rate, 
and therefore to the neutrino flux. As a consequence, a 
limit on l//e can be translated into a limit on the cosmic 
ray production in specific models, and may even challenge 
the paradigm that GRBs are the sources of the UHECR; 
see, e.g., Ref. [5]. It can be read off from Fig. |3] that 
IC59+40 already constrains the baryonic loading to be- 
low about 20. Future IceCube data may reach l//e ~ 2, 
which, in turn, means that about one third of the energy 
of the GRB goes into photons. If this fraction is unac- 
ceptably large for theoretical GRB models, they will be 
finally challenged. 

There are, however, some caveats in this analysis: For 
most of the bursts used for the IC40 analysis the redshift 
is actually not measured. In these cases, for long bursts 
in Ref. z = 2.15 has been used for the computation 
of the breaks, and L]^ = 10^^ erg s^^ for the compu- 
tation of /tt in Eq. (fTl). However, in a self-consistent 
numerical approach, it is clear that these two quanti- 
ties are correlated (c/., Eq. ^ for measured, with 
Li^° ~ (1 -I- z)E'^°/Tgo). We have in Fig. [2] followed 
the IceCube logic in these cases: we have computed the 
proton density from directly, and the photon density 
from the chosen (as it enters - using the cho- 
sen value of z in both cases. Note, however, that fixing 
z and computing the luminosity from F^ with Eq. ([2]) 
in a self-consistent approach leads to a prediction which 
is strongly dependent on the chosen "standard value" of 
z. In this case, choosing z ~ 2 clearly overestimates the 
neutrino flux, where we find that one burst with a large 
photon fluence dominates. On the other hand, the peak 
contribution to the diffuse flux may rather come from 
2: ~ 1 [15] . Because of this uncertainty, a conservative 
lower limit on the neutrino flux can be only obtained 
from the bursts with z measured, c/., lower dashed curve 
in Fig. [3j A similar logic could be applied to ty and F. 

In summary, we have revised the GRB fireball neutrino 
flux calculation from a particle physics perspective. Com- 
pared to earlier analytical computations, our numerical 
simulation takes into account the full spectral (energy) 
dependencies of the proton and photon spectra, as well 
as the cooling of the secondaries, flavor mixing, and addi- 
tional multi-pion, kaon, and neutron production channels 
for the neutrinos. We have found a significant devia- 
tion in the normalization of the neutrino flux prediction 
of about one order of magnitude, with a very different 
spectral shape peaking at slightly higher energies. We 
have shown from the re-computation of the analytical 
models where the discrepancies come from, and we have 
demonstrated that they all add up in the same direction. 
Note that the exact size of the corrections depends on the 
individual burst parameters, which we have taken into 



account. We have also demonstrated, by the reproduc- 
tion of the IC40 analysis, that our prediction is signifi- 
cantly below the current IceCube limit, which means that 
the conventional GRB fireball phenomenology is not yet 
challenged. Finally, we have quantified additional astro- 
physical and method-dependent systematical uncertain- 
ties in the computation for this particular model in order 
to determine a lower bound for the predicted neutrino 
fiux. We have demonstrated that the baryonic loading 
of the fireballs can be constrained with the full-scale ex- 
periment after about ten years to a level which will ex- 
ert significant pressure on the parameter space for GRBs 
as sources of the UHECR. Note, however, that we have 
only considered the simplest possible one-zone internal 
shock model, whereas the spectral shape may change in 
multi-zone models (see, e.g., Ref. [5D]) or taking into ac- 
count additional effects, such as the acceleration of the 
secondaries |21j . In addition, photospheric emission and 
magnetic reconnection models may lead to larger emis- 
sion radii, see, e.g., Ref. [22], and UHECR acceleration 
may also be possible at external shocks [531 121] ~ with 
potentially smaller neutrino fluxes. 
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